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Abstract 



We study macroscopic entanglement of various pure states of a one-dimensional iV-spin system 
with N 3> 1. Here, a quantum state is said to be macroscopically entangled if it is a superposition 
of macroscopically distinct states. To judge whether such superposition is hidden in a general state, 
we use an essentially unique index p: A pure state is macroscopically entangled if p = 2, whereas it 
may be entangled but not macroscopically if p < 2. This index is directly related to fundamental 
stabilities of many-body states. We calculate the index p for various states in which magnons are 
excited with various densities and wavenumbers. We find macroscopically entangled states (p = 2) 
as well as states with p = 1. The former states are unstable in the sense that they are unstable 
against some local measurements. On the other hand, the latter states are stable in the senses 
that they are stable against any local measurements and that their decoherence rates never exceed 
O(N) in any weak classical noises. For comparison, we also calculate the von Neumann entropy 
Sn/2{N) of a subsystem composed of N/2 spins as a measure of bipartite entanglement. We find 
that Sn/2(N) of some states with p = 1 is of the same order of magnitude as the maximum value 
N/2. On the other hand, S^-m^N) of the macroscopically entangled states with p = 2 is as small 
as O(logiV) <C N/2. Therefore, larger S^/ 2 (N) does not mean more instability. We also point 
out that these results are partly analogous to those for interacting many bosons. Furthermore, the 
origin of the huge entanglement, as measured either by p or S^/2(N), is discussed to be due to 
spatial propagation of magnons. 

PACS numbers: 03.67.Mn, 75.45. +j, 75.10.Jm, 03.65.Yz 
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I. INTRODUCTION AND SUMMARY 



Many-partite entanglement, i.e., entanglement in a system that is composed of many sites 



or par ties, has been attracting much attention recently 
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22j. It is known that the number of possible measures of 
entanglement grows dramatically as the number of sites is increased j^J. Different measures 
are related to different physical properties. Therefore, one must specify physical properties 
of interest in order to determine a proper measure or index. 

In this paper, we study macroscopic entanglement of various states in a quantum many- 
spin system. Here, a quantum state is said to be macroscopically entangled if it is a super- 
position of macroscopically distinct states (see Sec. 1111 A|) . Although such superposition is 
trivially recognized for some states (such as the 'cat' state), it is hard to find such superpo- 
sition by intuition for general states. In order to judge whether such superposition is hidden 
in a general state, we use an essentially unique index p, defined by Eq. (fTBjl. A pure state is 
macroscopically entangled if p — 2, whereas it may be entangled but not macroscopically if 
p < 2. Unlike many other measures or indices of entanglement, there is an efficient method 
of computing p for any given states [3j. 

It was shown by Shimizu and Miyadera (hereafter refereed as SM) that this index is 
directly related to fundamental stabilities of many-body states, i.e., to fragility in noises or 
environments and to stability against local measurements. That is, a state with p = 1 is not 
particularly unstable against noises in the sense that its decoherence rate does not exceed 
O(N) in any noises or environments, whereas the decoherence rate of a state with p = 2 
can become as large as 0(N 2 ) j^j. Furthermore, a quantum state with p = 2 is unstable 
against local measurements, whereas a homogeneous state with p — 1 is stable. 

We consider a one-dimensional A-spin system with A r >l, and calculate the index p for 
various pure states in which magnons are excited with various densities and wavenumbers. 
We find macroscopically entangled states (p = 2) as well as 'normal' states with p = 1 which 
are entangled but not macroscopically. According to SM, they are unstable and stable 
many-body states, respectively. 

For comparison, we also calculate the von Neumann entropy Sn/2{N) of a subsystem 
composed of A/2 spins as a measure of bipartite entanglement. We find that some states 
with Sn/2(N) = O(N), which is of the same order of magnitude as the maximumvalue N/2, 
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are 'normal' states in the sense that p = 1. On the other hand, some of other states, which 
are macroscopically entangled (p — 2), have much smaller value of Sn/2(N) of O(logiV). 

These results demonstrate that the degrees of entanglement are totally different if different 
measures or indices are used. Furthermore, stabilities of quantum states are not simply 
related to the degrees of entanglement: Different stabilities are related to different measures 
or indices. In particular, fragility in noises and the stability against local measurements are 
directly related to p, hence are basically independent of Sn/2(N). 

The results also demonstrate that states with huge entanglement, as measured by either 
p or Sn/2(N), can be easily constructed by simply applying creation operators of magnons 
to a ferromagnetic state, which is a separable state. Neither randomness nor elaborate tun- 
ing are necessary to construct states with huge entanglement from a separable state. This 
should be common to most quantum systems: By exciting a macroscopic number of elemen- 
tary excitations, one can easily construct states with huge entanglement. To generate such 
states experimentally, however, one must also consider the fundamental stabilities mentioned 
above: States with p = 2 would be quite hard to generate experimentally, whereas states 
with large Sn/ 2 {N) would be able to be generated rather easily. 

The present paper is organized as follows: In Sec.|nJ we shortly review physics of magnons, 
and present state vectors of many-magnon states under consideration. We explain the 
index p for the macroscopic entanglement, and present an efficient method of computing 
p in Sec. IIHI In Sec. IIV| we study macroscopic entanglement of many-magnon states by 
evaluating p. We study their bipartite entanglement in Sec. El for comparison purposes. 
Stabilities of the many-magnon states are discussed in Sec. IVII In Sec. IVH A| we point 
out that our results are analogous to those for interacting many bosons. Furthermore, we 
discuss the origin of the huge entanglement in Sec. IVII Bl 

II. MANY-MAGNON STATES 

In this section, we briefly review the physics of magnons in order to establish notations. 

A magnon is an elementary excitation of magnetic materials. It is a quantum of a spin 
wave that is a collective motion of the order parameter, which is the magnetization M. for 
a ferromagnet. 

For example, consider a one-dimensional Heisenberg ferromagnet which consists of N 
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spin-i atoms. Under the periodic boundary condition, its Hamiltonian is given by 



N 



H = -J^<?(Z) •<?(/ + !). 



(1) 



1=1 



Here, J is a positive constant, and <?(/) = (d x (l),a y {l),a z {l)), where a x (l),a y (l),a z (l) are 
Pauli operators on site I. We denote eigenvectors of a z corresponding to eigenvalues +1 and 
-1 by | |) and | j), respectively. One of the ground states of the Hamiltonian is | |® ), in 
which Ai points to the — z direction. The state in which one magnon with wavenumber k 
is excited on this ground state is 

N 

(2) 



1 - 

1=1 



where fr+(l) = (^(Z) + ia y (l))/2. The excitation energy of iV'fejjv) is easily calculated as 



E, 



k:N 



Usin 2 -. 



(3) 



It goes to zero as k — > because a magnon is a Nambu-Goldstone boson. The dispersion 
relation for small k is nonlinear because a magnon is a non-relativistic excitation. Because 
of the periodic boundary condition, k takes discrete values in the first Brillouin zone, — n < 
k < 7r: 



, 2vr . N N 

k = jrj U : integer, -— <]<—). 

It is conventional to define the 'creation operator' of a magnon with wavenumber k by 



(4) 



Ml 



The commutation relations are calculated as 



1 - 



e lkl a + (l). 



(5) 



Ml Ml 



M k ,Ml 



M k , M k < | = 0, 



-V, 

N ^ 



,i(k'-k)l 



(6) 
(7) 



i=i 



When the number m of magnons is much smaller than N, Eq. (0) can be approximated as 



N 



(8) 



Therefore, magnons behave as bosons when m <^ N. 

Using the creation operators, we define the m-magnon state with wavenumbers 
k 1: k 2 ,...,k m by 



m 



V l<>a-i<>b- 



i=l 



Here, l/^/nJn^. • ■ ■ is the usual normalization factor for bosons, where n u {y = a, b, ■ ■ • ) 
denotes the number of k^s having equal values, and Gk lt ...,k m -,N is an extra normalization 
factor which comes from the fact that magnons are not strictly bosons. Without loss of 
generality, we henceforth assume that 

h<k 2 < ...<k m . (10) 

When m <C N ', the magnons behave as bosons so that Gk u ...,k m ;N = 1 and 

{^k 1 ,k 2 ,-,k m ;N\^k' 1 ,k' 2 ,...,k' m ,N) = <Wi<W^, ' ' ' h m ,k' m (H) 

to a good approximation. On the other hand, the deviations from these relations become 
significant when m = O(N). 

An m-magnon state with a small density (m/N <^ 1) of magnons is an approximate en- 
ergy eigenstate. Although an m-magnon state with a large number (m = O(N)) of magnons 
is not generally a good approximation to energy eigenstate, such a state is frequently used 
in discussions on a macroscopic order because many magnetic phase transitions can be re- 

— * 

garded as condensation of O(N) magnons. For example, the state in which M. points to a 
direction with the direction vector (sin 9 cos a, sin 9 sin a, cos 9) can be described as 

\(9ar N ) = ^-cos-|T) + sin-||)J (12) 

N 

= ^2e- ima y/B^\^ {k=0)m . N ), (13) 

m=0 

where |'0(fe=o) r ";w) is the m-magnon state with k\ — . . . — k m — 0, and B m is the binomial 
coefficient; 

M:)K)>iP 

When 9 ^ n, B rn has a peak at m = iVcos 2 | = O(N), and thus a macroscopic number of 
magnons are 'condensed'. 

6 



Note that | 1® N ) and \ (8a)® N ) belong to the same Hilbert space because we assume that 
N is large but finite, although they will belong to different Hilbert spaces if we let N — > oo. 
For the same reason, all |^jfc 1 ,jfc 2 ,...,fc ro ;jv)' s belong to the same Hilbert space even if m = O(N). 



III. THE INDEX OF MACROSCOPIC ENTANGLEMENT 

In this section, we present the index of macroscopic entanglement, and an efficient method 
of computing it. We also explain its physical meanings by giving a few examples. Relation 
between this index and stabilities of many-body states will be explained in Sec. I VII 



A. The index p 



We are most interested in superposition of macroscopical 



been attracting much attention for many years 
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y distinct states, which has 



27, 28 . We note that such 



superposition was defined rather ambiguously. For example, the 'disconnectivity' defined 
in Ref. j^] is not invariant under changes of canonical variables, such as from the pairs of 
positions and momenta to the pair of a field and its canonical conjugate. Furthermore, in 
much literature the energy scale is not specified to determine the degrees of freedom involved 
in the superposition. However, the degrees of freedom, hence the disconnectivity, usually 
grows (decreases) quickly with increasing (decreasing) the energy scale under consideration 
|29[. On the other hand, SM proposed a new definition that is free from these ambiguities. 
We therefore follow SM. 

We first fix the energy range under consideration. For definiteness we assume that in 
that energy range the system can be regarded as N spin-| atoms, which are arranged on 
a one-dimensional lattice. We note that two states are 'macroscopically distinct' iff some 
macroscopic variable(s) takes distinct values for those states. As a macroscopic variable, 
it is natural to consider the sum or average of local observables over a macroscopic region 
[sflj ]. Since the average can be directly obtained from the sum, we only consider the sum 
in the following. That is, we consider additive operators [3l| . which take the following 
form: A = Y^Li^)- Here, a(l) is a local operator on site /, which, for the spin system 
under consideration, is a linear combination of the Pauli operators a x (l), cf v (1), a z (l) and the 
identity operator 1(1) on site I. Since we will consider all possible additive operators, we do 
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not assume that a(V) (/' 7^ I) is a spatial translation of a(l). 

Two states, and ^2), are macroscopically distinct iff some additive observable(s) A 
takes 'macroscopically distinct values' for those states in the sense that 

(Vi|i|V>i> - = O(N). (15) 

Therefore, if a pure state has fluctuation of this order of magnitude, i.e., if 6 A = 
[ (V; I Ai f Ai !?/>)] 1/2 = O(N) for some additive observable (s) A, where A A = A — (il)\A\ip), 
then the state is a superposition of macroscopically distinct states. On the other hand, if 



5 A = o(N) 23] for every additive observable A the state has 'macroscopically definite values' 
'or all additive observables. A typical magnitude of 5 A for such a state is 5 A = 0(N 1 ^ 2 ) 
. To express these ideas in a simple form, we define an index p for an arbitrary pure 
state \ip) by the asymptotic behavior (for large N) of fluctuation of the additive observable 
that exhibits the largest fluctuation for that state j^j]: 

sup^lAA+AilV') = 0{N P ). (16) 
AeA 

Here, A is the set of all additive operators. According to the above argument, \ip) is a super- 
position of macroscopically distinct states iff p = 2, and for pure states p is the essentially 
unique index that characterizes such a superposition. We therefore say that a pure state is 
macroscopically entangled iff p = 2. 



B. An efficient method of computing p 



It is well-known that many entanglement measures which are defined as an extremum 
are intractable Q S Q fl. In is an efficient -nethod of e 0mp ntn lg the 

index p [2]. We here explain the method briefly assuming an N spin-| system. 

Any local operator a(l) on site / can be expressed as a linear combination of 
cr x (l), d~ y (l), a z (l) and 1(Z). Since the identity operator 1 does not have fluctuation for any 
state, we can limit ourselves to local operators that are linear combinations of Pauli opera- 
tors. Therefore, an additive observable in question generally takes the following form; 

N N 



.4 



1=1 



1=1 a=x,y,z 
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where complex coefficients (see Sec. IIV Ej) . Since the local operators should not 

depend on iV (because otherwise A would not become additive), c^'s should not depend 
on N, hence the sum ^ J2 a \ c ai\ 2 is 0(N). Since we are interested only in the power p of 
AA^AAlifj) = 0(N P ), we can normalize c a i without loss of generality as 

N 

E E m 2 = n - (is) 

1=1 a=x,y,z 

The fluctuation of A for a given state \ip) is expressed as 

(V|AitAA|V) = £ c^K^,, (19) 

where V a i,pv is the hermitian matrix defined by 

= (V|Aa- a (Z)Ad- /9 (Z')|V>, (20) 

which we call the variance- covariance matrix (VCM) for It is seen from Eq. f!19|) that 
eigenvalues of this matrix are non- negative, and that (ip\ AA^ AA\ip) takes the maximum 
value when c a i is an eigenvector of the VCM corresponding to the maximum eigenvalue 
e ma x- By taking c a i of Eq. (flT?j) as such an eigenvector, we obtain 

sup^|AitAi|^) =e max iV. (21) 
AeA 

Therefore, e max is related to p as 

e max = 0(iVP- 1 ). (22) 

For example, p = 1 if e max = 0(1) whereas p = 2 if e max = 0(N). 

Note that we can evaluate p of a given state using this method in a polynomial time of 
the number iV of spins, because we have only to calculate the maximum eigenvalue of a 
VCM, which is a 3N x 3N matrix. 

C. Examples of macroscopically entangled states 

The iV-spin GHZ state, or the 'cat' state, |GHZ) = (| [® N ) + \ | 07V )) , violates a 
generalized Bell's inequality by a macroscopic factor j^. The index p correctly indicates that 
this state is macroscopically entangled, p = 2 2]. In contrast, Sn/ 2 (N) (which is defined by 
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Eq. (jHUj) below) of this state is extremely small; Sat/2 (AQ — 1. It may be intuitively trivial 
that this state is macroscopically entangled. However, intuition is useless for more general 
states such as the following examples. The greatest advantage of using p is that it correctly 
judges the presence or absence of macroscopic entanglement for any complicated pure states. 

For example, it was recently shown [2] that the quantum state of many qubits in a 
quantum computer performing Shor's factoring algorithm is transformed in such a way that p 
is increased as the computation proceeds, and the state just after the modular exponentiation 
processes, 

2*1 -1 

| ME) = -= l a )il x ° mod M >2 ? (23) 
V2 Nl a=Q 

is a macroscopically entangled state. Here, | • ■ ■ )i (| • • ^2) represents a state of the first 
(second) register, Aq (21ogM < N% < 21ogM + 1) denotes the number of qubits in the first 



register, x is a randomly taken integer, and M is a large integer to 



:>e factored. This state 
Although presence 



2 ] 



was shown to play an essential role in Shor's factoring algorithm 
of entanglement in this state is obvious, the presence of macroscopic entanglement was not 
revealed until an additive operator whose fluctuation is 0(N 2 ) was found in Ref. Q]. 

Another example is entanglement of ground states of antiferromagnets, which has recently 
been studied by many authors using various measures (f| Q, 0, 0, Ell- It is well-known 
that the exact ground state |Gaf) of the Heisenberg antiferromagnet on a two-dimensional 
square lattice of a finite size is not the Neel state but the symmetric state that possesses 



all the symmetries of the Hamiltonian 



3- 



We here point out that I Ga f) is entangled 



macroscopically, p = 2. In fact, the ground state has a long-range order jsTT ] . 

(GafKM^IGaf) ~ 0.117Af 2 + 1.02ATi , (24) 

where = Yl^Li(~ is the staggered magnetization (a = x,y,z). On the other 
hand, (Gaf|M?|Gaf) = by symmetry. Therefore, the order parameter of the antifer- 
romagnetic phase transition exhibits a huge fluctuation, 

(GafKAM^IGaf) = 0(N 2 ). (25) 



This shows that p = 2 for |Gaf)- Note that such a macroscopically entangled ground state 
appears generally in a finite system that will exhibit ajphase transition as A^ — > 00 if the order 
parameter does not commute with the Hamiltonian 0,0]. For example, the ground state of 
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interacting bosons |39l . |40( , for which the order parameter is the field operator of the bosons, 
is entangled macroscopically. Moreover, the ground state of the transverse Isin g m odel, 
whose entanglement has recently been studied using various measures J^, Q> 12 J3, Q , 
also has p = 2 when the transverse magnetic field is below the critical point. 

As demonstrated by these examples, the index p captures the presence or absence of 
certain anomalous features, which are sometimes hard to find intuitively, of pure quantum 
states in finite macroscopic systems. Furthermore, as we will explain in Sec. IVI1 p is directly 
related to fundamental stabilities of many-body states. 



IV. MACROSCOPIC ENTANGLEMENT OF m-MAGNON STATES 

We now study macroscopic entanglement of magnon states with various densities and 
wavenumbers by evaluating the index p. 



A. States to be studied 



Most relevant parameters characterizing the magnon states are the number of magnons, 
m, and the wavenumbers of magnons. Because of the Z 2 symmetry, we assume that 1 < 
m < N/2 without loss of generality. Furthermore, we assume that N is even in order to 
avoid uninteresting complications. 

Since we are interested in the asymptotic behavior for large N, only the order of mag- 
nitudes of these parameters is important. We therefore consider the following three cases 
41|: 

(a) m = 0(1). 

(b) m = 0(N) and all magnons have different wavenumbers from each other, continuously 
occupying a part of the first Brillouin zone. Because of the translational invariance 
of the system in the fc-space, it is sufficient to calculate the case where the magnons 
continuously occupy the first Brillouin zone from the bottom, i.e., their wavenumbers 
are 0, • • • , respectively. 

(c) m = O(N) and all magnons have equal wavenumbers k. Because of the translational 
invariance of the system in the fc-space, we can take k = without loss of generality. 
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Furthermore, a small number (= 0(1)) of magnons with arbitrary wavenumbers may be 
added to these states. It is expected and will be confirmed in the following that the addition 
does not alter the value of p. 

B. Case (a) 

In Fig. Q we plot numerical results for e max of two-magnon states as functions of N. The 
result for k\ = k<i can also be obtained from the analytic expression, Eq. (J37J). These results 
show that excitation of a small number (0(1)) of magnons on the ferromagnetic ground 
state | 1® N ), which is a separable state, does not change the value of p. It is thus concluded 
that magnon states for case (a) are not macroscopically entangled. 

7 
6 
5 

x 4 

E 

° 3 
2 
1 


10 20 30 40 50 60 

N 

FIG. 1: The maximum eigenvalue e max of the VCM of two-magnon states with wavenumbers k\ 
and k% as functions of the number N of spins. Because of the translational invariance of the system 
in the fc-space, we take k\ = without loss of generality. The solid line represents the analytic 
expression e max = 1 + (5iV — 12) /N, Eq. (|27|). which assumes that all wavenumbers are equal. 

C. Case (b) 

To investigate p for case (b), we evaluate e max for various magnon densities assuming that 
all magnons have different wavenumbers from each other, continuously occupying the first 
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Brillouin zone from the bottom. The results are plotted in Fig. El for m = N/2, N/4, and 
N/6 as functions of N. It is seen that e max ~ constant, hence p — 1. We also confirmed (not 
shown in the figure) that addition of small number of magnons with arbitrary wavenumbers 
does not alter the value of p. We thus conclude that magnon states for case (b) are not 
macroscopically entangled. 
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FIG. 2: The maximum eigenvalue e max of the VCM of m-magnon states with m = N/2, N/4, and 
iV/6 as functions of the number iV of spins. The wavenumbers of magnons are all different taking 
the values 0, ±2n/N, ±4tt/N, . . ., respectively, i.e., the first Brillouin zone is continuously occupied 
from the bottom. 



D. Case (c) 

If the wavenumbers of all magnons are equal, we can calculate e max analytically as fol- 
lows. Since we can take k = by symmetry, we calculate the VCM of |V , (fc=o) m ; v)- From 
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calculations described in Appendix A, we obtain the VCM and the maximum eigenvalue as 



V, 



j -j 



1 


(j = f 


1 < j < 2N) 


i-wi 


(j = f 


2N + 1 < j < 3N) 


Wi 


(j ^ f 


l<j,j'<N) 


Wi 


(j + f 


N + l<j,j'<2N) 


w 2 - wi 


(j ^ f 


2N + l<j,j'<3N) 


-iW 3 


(j = f 


— N, N + l<f< 2N) 


iW 3 


(j = f 


+ N, l<j'< N) 





(others 


), 



(26) 



1 + (N - 1)W 1 + W 3 

2mN - 2m 2 + N - 2m 
1 + N ■ 



(27) 



where W\, W2, and W 3 are defined by Eqs. (|A2j) . ()A3|) . and (JAJJ), respectively. We therefore 
find that e max = 0(N) for m = 0(N), hence p = 2. 

The solid line in Fig. ^represents the analytic expression for e max , Eq. (|27|h for m = N/2. 
We also plot numerical results for the cases where the wavenumbers of one or two magnons 
are different. It is seen that e max becomes smaller in the latter cases, as we have seen a similar 
tendency in Fig. [TJ However, e max = 0(N) and thus p = 2 in all three cases in Fig. El We 
therefore conclude that magnon states for case (c) are macroscopically entangled. 



E. Additive operator with the maximum fluctuation 

For a given state \ip), we can obtain the additive operator A max that has the maximum 
fluctuation (5A max ) 2 = {ip\ AA max AA max |-^) = Ae max for that state by inserting the eigen- 
vector of the VCM belonging to the maximum eigenvalue e max into Eq. (JT7J). However, A max 
is generally non-hermitian because the eigenvector is generally complex. A non-hermitian 
operator A can be decomposed into the sum of two hermitian operators A 1 and A" as 
A = A' + iA". If A' and A" commute with each other, they can be measured simultaneously 
with vanishing errors. Since the values of A have one to one correspondence to the pairs 
of the values of A' and A", one can measure A by simultaneously measuring A' and A" if 
[A', A"] = 0. Note that in macroscopic systems [A', A"} ~ to a good approximation for 
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FIG. 3: The maximum eigenvalue e max of the VCM of iV/2-magnon states as functions of the 
number N of spins. Most magnons have equal wavenumbers, i.e, most magnons are 'condensed.' 
The solid line represents the analytic expression c max — 1 ~t~ Ecj. (|27|) . which assumes that all 
wavenumbers are equal. The circle and crosse represent numerical results for the cases where the 
wavenumbers of one or two magnons are different. Because of the translational invariance in the 
k-sp&ce, we take the wavenumber of the condensed magnons as without loss of generality. 

any additive operators A' and A" because [(A'/N), (A"/N)] is at most 0(1/N) ~ 0. There- 
fore, in macroscopic systems non-hermitian additive operators can be measured to a good 
accuracy. Hence, A max can be measured even if it is non-hermitian. One can also construct 
the hermitian additive operators A' max = (i max + and A'^ = (A max - i^ ax )/2i, 

which are the 'real' and 'imaginary' parts, respectively, of A max . Using the triangle inequal- 
ity ||AA nax |^)|| < ||Av4^ ax |^)|| + ||AA^ ax |^)||, we can easily show that either (or both) of 
SA' mSuX or SA'^ ax is of the same order as <5A max . 

For \ip(k=o) m ;N) with m = O(N), for example, the eigenvector belonging to the maximum 
eigenvalue (|27j) is 

N N N 

(28) 




which gives the maximally fluctuating additive operator as 

1 N 



N 



(29) 



i=i 



i=i 



15 



for which (SA maiX ) 2 = 0(N 2 ). Although this operator is not hermitian, it can be measured to 
a good accuracy if N 3> 1. Or, let us define hermitian operators A' max = (A max + v4 max )/2 = 

|Z)£=i^(0 and ^max = (^max - ^Lax)/ 2 * = | E/Ii <^ (0 • Since 1^=0)-;^) is symmetric 
under rotations about the z axis, we can show that {5A' max ) 2 = (<W max ) 2 = 0(N 2 ) in this 
case. It is worth mentioning that A max corresponds to the eigenvector belonging to the 
second largest eigenvalue e 4 , which is given by Eq. (|A4|) and is of 0(N). 

V. BIPARTITE ENTANGLEMENT OF m-MAGNON STATES 

For a comparison purpose, we now calculate the degree of bipartite entanglement of 
magnon states that have been studied in the previous section. For a measure of bipartite 
entanglement, we use the von Neumann entropy of the reduced density operator of a sub- 
system. That is, we halve the iV-spin system and evaluate the reduced density operator 
Pn/2{N) of one of the halves. The von Neumann entropy is defined by 

S N/2 (N) = -Tr [p N/2 (N) \og 2 p N/2 (N)] , (30) 

which ranges from to N/2. Although Sn/2(N) for the case where all wavenumbers are 
equal was discussed by Stockton et al. we here evaluate Sn/ 2 (N) systematically for all 
the three cases listed in Sec. HVH 

A. Case (a) 

To evaluate Sn/2{N), we halve the iV-spin system into two subsystems A and B. Ac- 
cordingly, we decompose \ipk 1 ,k 2 ,--;k m ;N) into the sum of products of \ipk 1 ,k 2 ,---;N/2)'s of A and 
B. 

When all wavenumbers are different from each other, an m-magnon state can be decom- 
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posed as 



i=i 

Gk lt k 2 ,...,k m ;N ( n -\ I |®jV/2\|„ /, \ 
^= ^fci.-.fcn.SiV/al i )Wkx,...,k m ;N/2) 



rn 



+ ^ ^N/iQ^Q^ ^ ... Mm . N/2 \^N/2) lK...A,..,^/2> 
i=l 

m— 1 m 

+ ejfcl ^ /2+ifcj7V/2G fci J ;Af/2 G A,. 1 1 ;... ,4... ,fc ^jA^l^i^/^ l^fe!,...,^,...,^,...,^;^) 



i=l j =i+l 



+ e ^^ + ...+^^ G -i . . ^I^,...,^)! I^ /2 > ) , (31) 



%,-,fc m ;iV/2 

where ~ denotes absence, and the prefactor l/y/2™ comes from the prefactor 1/y/N in 
Eq. (jSJ). When m = 2, for example, 

+e iklN/2 G~^. N/2 G~^ N/2 \^ kl . N / 2 ) \ipk 2 ;N/2) + e lk2N/2 G~^. N/2 G~^ N/2 \^ k2 . N/2 )\tlj kl . N / 2 } 
+ ^ lJV/2+ifc2jV/2 ^ fc2;2V/2 |^ llte;J v/2)| i^ /2 >) , (32) 

which means that the state is a superposition of the following four (= 2 2 ) states: (i) both 
magnons are in subsystem B, (ii) the magnon with k\ is in A whereas the magnon with &2 
is in B, (iii) the magnon with k 2 is in A whereas the magnon with k\ is in B, and (iv) both 
magnons are in A. As N is increased in decomposition (J31|) (while m is fixed), all G"s — > 1 
and 2 m vectors on the right-hand side tend to become orthonormalized. This means that 
decomposition (J3~T|) becomes the Schmidt decomposition, in which the Schmidt rank is 2 m 
and all the Schmidt coefficients are equal (except for the phase factors). We thus obtain 

2™ / -i \ 2 / -, \ 2 

= m, (33) 



Ito Sk/2 (JV) = -£U=) log 2 (-U 
■ WH^ i=i vvz y vvz y 



(m: fixed) 

i.e., Sn/2{N) = 0(1) for fixed m. Note that m is the maximum value of Sn/2{N) among 
states whose Schmidt rank is 2 m . 

When some of the wavenumbers are equal, the Schmidt rank becomes smaller because 
magnons having equal wavenumbers are indistinguishable. For example, if k\ = k 2 = k, the 
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two-magnon state 

hfc,*;*> = ^{Mlf\ 1® N ) (34) 

is decomposed as 

;JV/2/ 

+ ^G-i ;7V/2 |^ /2 )||^ /2 )). (35) 

In contrast to Eq. (|32p. this is the superposition of the following i/iree (< 2 2 ) states; (i) both 
magnons are in 5, (ii) one magnon is in A and the other is in B, and (iii) both magnons 
are in A. The Schmidt rank is thus decreased. Furthermore, the Schmidt coefficients do not 
take the same value. As a result of these, Sn/2(N) becomes smaller than Eq. (J33J); 

lim S N/2 (N) < m, (36) 
(m: fixed) 

from which we again have Sn/2{N) = 0(1). 

Therefore, we conclude that the bipartite entanglement of magnon states in case (a) is 
small in the sense that 

S N/2 (N) = 0(1). (37) 

As a demonstration, we plot numerical results for Sn/2{N) as functions of N in Fig.|U for 
three-magnon states in the following three cases; (i) three magnons have different wavenum- 
bers fci, &2, Afe, (ii) two magnons have equal wavenumbers k\ = k 2 whereas one magnon has 
another wavenumber k%, and (iii) three magnons have equal wavenumbers k± = k 2 = ks- For- 
mulas (|3*3*|) and (J3fi|) are confirmed. Furthermore, it is seen that the departure of magnons 
from ideal bosons becomes significant for small N, and that Sn/2{N) approaches the lim- 
iting values for iV — > oo from below. This may be understood from the discussions of the 
following subsections. 

Note that the result ()37|) agrees in some sense with the result of Sec. II V Bl in which we 
have seen that the states are not macroscopically entangled. However, we will see in the 
following that such a simple agreement is not obtained in cases (b) and (c). 
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FIG. 4: The von Neumann entropy of a subsystem, S^/ 2 (N), of three-magnon states with wavenum- 
bers k\, k 2 , and /C3 as functions of the number N of spins. Because of the translational invariance 
of the system in the /c-space, we take k 2 = without loss of generality. 

B. Case (b) 

In case (b), the previous argument on limjv^oo Sn/2(N) does not hold because the de- 
parture of magnons from ideal bosons is significant when m = O(N). In fact, the vec- 
tors in Eq. ()31|) do not become orthonormalized as iV — > oo. Hence, Eq. (|31|) does not 
become the Schmidt decomposition, and it can be further arranged until it becomes the 
Schmidt decomposition. Therefore, we expect that the Schmidt rank is less than 2 m , and 
limjv^oo S N/2 (N) < m. 

To see more details, we have calculated Sn/2(N) numerically. The results are plotted as 
functions of N in Fig. El for case (b) with m = N/2,N/4, and N/6. It is found that the 
results are well approximated by the straight lines, 



S N/2 (N) = aN + b, 



(38) 



which are also displayed in Fig. |SJ The parameters a and b are determined by the least 
squares, whose values are tabulated in Table HJ Since < a < m/N, we find that Sn/ 2 (N) 
is less than, but of the same order of magnitude as, the maximum value N/2; 



S N/2 (N) = 0(N). 
19 



(39) 



TABLE I: The parameters a and b, which are calculated with the least squares, of the regression 
line Eq. Q38[) for m-magnon states of Fig. [21 



m 


a 


asymptotic standard error 


b 


asymptotic standard error 


N/2 


0.36 


± 0.009 


0.27 


± 0.113 


N/4 


0.21 


± 0.002 


0.29 


± 0.027 


N/6 


0.15 


± 0.003 


0.15 


± 0.044 



We thus conclude that the bipartite entanglement of magnon states in case (b) is extremely 
large. This should be contrasted with the result of Sec. IIV CI according to which these states 
are not macroscopically entangled. 
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FIG. 5: The von Neumann entropy of a subsystem, S^/z^N), of m-magnon states with m = 
N/2,N/A, and iV/6 as functions of the number iV of spins. The wavenumbers of magnons are 
all different taking the values 0, ±2n/N, zt&ir/N, . . ., respectively, i.e., the first Brillouin zone is 
continuously occupied from the bottom. The lines represent the regression lines calculated with 
the least squares. 
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C. Case (c) 



We finally consider Sn/2{N) in case (c). When all wavenumbers are equal to zero, the 
m-magnon state W>(k=o) m ;N) becomes identical to the "Dicke state" that was discussed by 
Stockton et al. [15]. According to their result, 

S N/2 (N) = 0(logN) (40) 

when m = O(N). Because of the translational invariance, this result also holds for \ip(k) m -,N) 
with other values of k. Since 1 < O(logiV) <C O(N), we find that SV/ 2 (iV) is slightly larger 
than that of case (a), but much smaller than that of case (b). We therefore conclude that 
the bipartite entanglement of magnon states for case (c) is small. This should be contrasted 
with the result of Sec. II V Dl according to which these states are macroscopically entangled. 

VI. STABILITIES AND ENTANGLEMENT 

It may be expected that a quantum state with larger entanglement would be more unsta- 
ble. This naive expectation is, however, quite ambiguous for many-body systems. First of 
all, the degree of entanglement depends drastically on the measure or index used to quantify 
the entanglement, as we have shown above. Furthermore, 'stability' can be defined in many 
different ways for many-body states. 

n 

SM considered the following two kinds of stabilities One is the stability against 
weak perturbations from noises or environments: A pure state is said to be fragile if its 
decoherence rate behaves as ~ KN 1+5 when perturbations from the noises or environments 
are weak, where 5 is a positive constant. Such a state is extremely unstable in the sense 
that its decoherence rate per spin increases as ~ KN S with increasing N, until it becomes 
extremely large for huge N however small is the coupling constant between the system and 
the noise or environment. SM showed that pure states with p = 1 never become fragile in 
any noises or environments, whereas pure states with p = 2 can become fragile, depending on 
the spectral intensities of the noise or environment variables. The other stability considered 
by SM is the stability against local measurements: A state is said to be stable against local 
measurements if an ideal (projective) measurement of any observable at a point / does not 
alter the result of measurement of any observable at a distant point /' for sufficiently large 
\l — l'\. SM showed that this stability is equivalent to the 'cluster property', which is closely 
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related to p, if the cluster property for finite systems is properly defined. For example, a 



42| state with 



state is unstable against local measurements if p = 2, whereas a homogeneous 
p = 1 is stable. 

We have shown that p = 1 for magnon states of case (b). Therefore, these states never be- 
come fragile in any noises or environments, and they are stable against local measurements. 
Since we have also shown that the bipartite entanglement of these states is extremely large, 
we find that the bipartite entanglement is basically independent of these fundamental stabil- 
ities. The same conclusion was obtained for chaotic quantum systems by two of the authors 
0|: They showed that p = 1 for almost all energy eigenstates of macroscopic chaotic systems 
whereas their bipartite entanglement is nearly maximum. 

We have also shown that p = 2 for magnon states of case (c). Therefore, these states can 
become fragile, depending on the spectral intensities of the noise or environment variables. 
Furthermore, these states are unstable against local measurements. Since we have seen 
that the bipartite entanglement of these states is small, we find again that the bipartite 
entanglement is basically independent of these fundamental stabilities. To understand the 
physics of this conclusion, the following simple example may be helpful. The W-state, 

iw) = 4fEi^ _1) ti 0M ). ( 41 ) 



has the same value of Sn/2{N) as the iV-spin GHZ state, i.e., Sn/2(N) = 1. On the other 
hand, p — 1 for |W) whereas p = 2 for |GHZ). As a result, the decoherence rate of |W) never 
exceeds O(N) by any weak classical noises, whereas the decoherence rate of |GHZ) becomes 
as large as 0(N 2 ) in a long-wavelength noise. Furthermore, |W) is stable against local 
measurements, whereas |GHZ) is unstable. These results are physically reasonable because 
the W-state is nothing but a one- magnon state (with k — 0), which can be generated easily 
by experiment: Such a state does not seem very unstable. (See also Sec. IVIIBl ) 

It should be mentioned that another stability was studied by Stockton et al. 
for a special state of case (c), i.e., for a Dicke state which in our notation is writ- 
ten as \?P(k=o) N / 2 ;N)- They showed that a bipartite entanglement measure of Pn-n> = 
Tr N i(\ip( k=0 jN/2. N )(ip( k= tyN/2. N \) decreases very slowly as N' increases. Here, Tr N i means 
'trace out N' spins'. They thus concluded that the state is robust. Although their conclu- 
sion might look contradictory to our conclusion, there is no contradiction. The stability 
(robustness) as discussed by Stockton et al. is totally different from the fundamental sta- 
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bilities that are discussed in the present paper. The Dicke state is 'robust' in the sense of 
Stockton et al, whereas in the senses of SM the state is 'fragile' in noises or environments 
and 'unstable' against local measurements. This demonstrates that stability can be defined 
in many different ways for many-body states. 



VII. DISCUSSIONS 



A. Relation to Bose-Einstein condensates 



The 'cluster property' which is closely related to the index p, of condensed states of 
interacting many bosons was previously studied in Ref . . It was shown there that p = 2 
for the ground state \N, G), which has a fixed number N of bosons, if N is large enough 
to give a finite density for a large volume. Since magnons are approximate bosons, magnon 
states of case (c) may be analogous to this state. Although deviations from ideal bosons 
become significant in case (c), the deviations may be partly regarded as effective interactions 
among magnons. This analogy intuitively explains our result that p = 2 for magnon states 
of case (c). 

It was also shown in Ref. ^Qj that p — 1 for a generalized coherent state |a,G), which 
was called there a coherent state of interacting bosons. This result may also be understood 
intuitively on the same analogy. That is, \a, G) may be analogous to the state of Eq. (J13j) . 
which has p = 1 because it is a separable state as seen from Eq. (|12jl. Therefore, by analogy, 
p should also be unity for \a, G), in consistency with the result of Ref. j^jj], although \a, G) 
is not separable. 

Analogy like these may be useful for further understanding of systems of interacting many 
bosons and of many magnons. 



B. What generates huge entanglement? 

We have shown that states with huge entanglement, as measured by either p or Sn/2{N), 
can be easily constructed by simply exciting many magnons on a separable state. We now 
discuss the physical origin of this fact. 

The most important point is that a magnon propagates spatially all over the magnet j^J . 
By the propagation, quantum coherence is established between spatially separated points 
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|44( | . Therefore, by exciting a macroscopic number of magnons, one can easily construct 
states with huge entanglement. 

Note that this should be common to most quantum systems, because Hamiltonians of 
most physical systems should have a term which causes spatial propagation. For example, 
such a term includes the nearest-neighbor interaction of spin systems, the kinetic-energy 
term of the Schrodinger equation of particles, the term composed of spatial derivative of 
a field operator in field theory. Therefore, excitation of a macroscopic number of elemen- 
tary excitations generates huge entanglement. Neither randomness nor elaborate tuning is 
necessary. 

This observation will be useful for theoretically constructing states with huge entangle- 
ment. Experimentally, on the other hand, the stability should also be taken into account 
because unstable states would be hard to generate experimentally. We thus consider that 
states with p = 2 should be much harder to generate experimentally than states with large 
Sn/2{N). In other words, a state with large Sn/2(N) would be able to be generated rather 
easily, e.g., by exciting many quasi-particles in a solid. In this respect, a naive expectation 
that states with large entanglement would be hard to generate experimentally is false: It 
depends on the measure or index that is used to quantify the entanglement. 
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APPENDIX A: CALCULATION OF THE VCM AND ITS EIGENVALUES 

The state vector of the m-magnon state with k\ — ■ ■ ■ — k m — can be written as 

|-0( fe =O)-;iV> = E E E • • • * + (!m)\l 9N ), 

^ N m h h{>h)h(>h) l m (>l m -l) 

where ArC m = . Since the VCM is hermitian, we have only to calculate the fol- 

\m J 

lowing correlations; (Aa x {l)Aa x (l')), (Aa x (l)Aa y (l')), (Aa x (l)Aa z (l')), (Aa y (l)Aa y (l')), 
(Aa y (l)Aa z (l')), {Aa z (l)Aa z (l')), where (•) stands for (V'^o)" 1 ;^! • IV'C^o)™;^)- Since 
IV'^owjv) is an eigenvector of exp(— id J2f=i °~z(l)), the state vector is invariant under a 
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rotation about z-axis. Therefore (Aa y {l)Aa y {l')) = {Aa x {l)Aa x {l')) and (Aa y {l)Aa z {l')) = 
(Aa x (l)Aa z (l')). Thus we calculate only 

(Aa x (l)Aa x (l')), (Aa x (l)Aa y (l')), (Aa x (l)Aa z (l')), (Aa z (l)Aa z (l')). 

We note that (<7 X (Z)) = (<r y (l)) = by symmetry, and 

1 N — 2m 
= -^riN-iCm-i - N-lC m ) = — - = -W 3 . (Al) 

When / = I', we easily obtain {a x {l)a x {l)) = (a z {l)a z {l)) = 1, ((T x (l)fr z (l)) = 0, and 

(a x (l)a y (l)) = i(a z (l)) = -iW 3 . 

Therefore, 

(Aa x (l)Aa x (l)) = 1, 
(Aa x (l)Aa y (l)) = -iW 3 , 
(Aa x (l)Aa z (l)) = 0, 
(Aa z (l)Aa z (l)) = 1-Wl 

When I 7^ /', we note that 

l V 

V NC m (ij( k= o) m . N \a x (l)d- x (l') \ . . . I ... t •••) = !, 

i v 

V NC m (4>( k =0)™;N\crx(l)c r x(l') \ ■ ■ ■ T ••■ I •••) = !> 

Z V 



V NC m (4)( k= Q)rn. N \a x (l)a x (l') \ . . . I ... I . . .) = 0, 

/ v 

V NC m (i)( k=0 )m; N \a x (l)cr x (l') \ . . . t ••• T • • •) = 0, 



where | . . . j ... j . . .) is a state vector in which m spins including Z'-th spin are up, 
whereas N — m spins including Z-th spin are down. We thus obtain 

/a - /,\a- n'\\ 2 N-2C m -i 2m(N-m) 

M|A "- (I|) = — = liM 5 »" (A2) 
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Furthermore, since 



V NC m {ljj(k=0)™;N\&x(l)(Ty(l')\ 

V NC m (ip^ k=0 )m. N \a x (l)a y (l')\ 



• I ••• T 



= i, 



• T 



V AfCm(V ; (fc=0) m ;Af| ( 5"a;(O^J/(^)l 
\/ NC m {ljj(k=0) m ;N 1 6x(0^S/(OI 



..) = o, 



. t ••• T •••) = o, 



we obtain 



(Aa x (l)Aa y (l>)) 



1 



(«Af-2C m -l — iN-2C m -i) — 0. 



It is obvious that (<3" a; (/)(3" z (/ / )) = because \ip(k=o) m ;N) is a linear combination of vectors 
whose m spins are up and N — m spins are down. Therefore 



(Aa x (l)Aa z (l')} = 0, 



Finally, since 



V AfCm('0(A:=O) m ;Ar|'5"z(O^(^)l 



V ivC m ('0(fe=O)'™;Ar|^(O^(OI 



i 



i ...) = 



V NC m (lp(k=0) m ;N\c l 'z(l)d'z(l') \ ■ ■ ■ T ••• T 



-1, 
-1, 
1, 
1. 



we obtain 



(a z (l)a z (l')) 



1 



2Ar_ 2 C m _i + N-2C m + ^-2^778-2) 



N 2 - 4mN - N + Am 2 



N(N-l) 



= W 2 



(A3) 



Therefore, 



{Aa z (l)Aa z (l')) = W 2 -W; 
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Combining these results, we obtain the VCM as Eq. (|26|) . or 
/ 



V 





Wi 


-iW 3 












1 





-iW 3 







iW 3 





1 


w x 










iW 3 


Wi 


1 

















i-wi 


w 2 - w 2 












w 2 - w 2 


i-wi 



\ 



We can calculate the eigenvalues 6j and the numbers Mj of the corresponding eigenvectors 
as 

Am(N — m) 

e 2 = 1 — W% + (N — 1)(W 2 - Wl) = 0; M 2 = 1, 



e 3 = 1 + W 3 + (N - l)Wa 



2N -2m + 2mN - 2m 2 

N ; 



M 3 



(A4) 



2m + 2mN - 2m 2 , , 
e 4 = 1 - W 3 + (iV - l)Wt = ; M 4 = 1, 

2N 2 -2N - 4mN + 2m + 2m 2 
e 5 = l + W 3 -W 1 = jv(iV-l) ; Af 5 = JV - 1 > 

e 6 = 1 - W 3 - W x = 2 ^~^y M e = N-l. 
The largest one is e 3 , which degenerates with when = 2m. We thus obtain Eq. (|27jl . 
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